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Abstract 

We assume that all quark and lepton 3x3 mass matrices which ap- 
pear in the standard model lagrangian (after spontaneous symmetry 
breaking) with neutrinos being treated as Dirac particles have the 
triangular form. Such matrices have not only less non-zero elements 
(tree of them are equal to zero) but also lead to very asymmetrical 
decomposition into one diagonal and two unitary matrices for quarks 
and leptons. We also assume that unitary matrices which transform 
flavor into definite mass states for right handed components (weakly 
non-interacting) in the same weak isospin doublet are equal. Using all 
available experimental data on quark and lepton masses and mixing 
angles, treating in the universal way quarks and leptons, we determine 
the triangular mass matrices for up and down type quarks, neutrinos 
and charged leptons and as a consequence mixing matrices for left- 
handed and right-handed components. As the result of the fit we get 
predictions for the neutrino masses including smallest neutrino mass. 
The calculations without CP violation and with inclusion of this effect 
in quark sector are also presented. 
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We have not only an information from the experiment about quark masses 
and Cabibbo-Kobayashi-Maskawa (CKM) fIJ mixing matrix for quarks j2] 
but now (after recent neutrino oscillation experiments [3]) one also knows two 
mass square differences for neutrinos and with less precision Maki-Nakagawa- 
Sakata (MNS) jl] mixing matrix for leptons [3J E]. Both mixing matrices 
(CKM and MNS) are products of matrices for left handed components for 
type u and type d quarks and for charged leptons and neutrinos. Knowing 
these matrices from experiment we do not know what is the connection be- 
tween left and right handed components of the flavor states and states with 
definite mass for up, down quarks, neutrinos and charged leptons. It was 
pointed out that this creates a problem when looking for patterns of masses 
and symmetries in neutrino and quark sectors |Hj • On the other hand if 
one knew from the lagrangian, after spontaneous symmetry breaking, mass 
matrices (they seem to be basic quantities) for quarks and leptons one could 
find by decomposition in two unitary matrices and diagonal matrix all the 
interesting mixing matrices. Then one is able to calculate CKM and MNS 
mixing matrices. 

The aim of our paper is, treating neutrinos as Dirac particles, to give the 
explicit example of possible mass matrices for u, d, v and l~ giving mixing for 
left handed and right handed components between states with definite flavor 
and definite mass. We will assume that the mass matrices for quarks and 
leptons gotten from standard model after spontaneous symmetry breaking 
are in form of triangular matrices. With additional assumptions discussed 
later we determine them using available experimental information on quark 
and lepton masses and CKM and MNS mixing matrices. We consider at the 
begining case without CP violation and then include this effect in the quark 
sector. 

Let us introduce the notation for quarks. The SU (3) * SU (2) * £7(1) gauge 
invariant Yukawa interactions for quarks are given by 

-Ly = Qi(Y u )ijURjH^ + Q t (Y d ) t:j d Rj H + h.c, (1) 

where Qj denote the SU (2) doublets of left-handed quarks and u R , d R are the 
right-handed up and down-type quarks respectively. The Yukawa couplings 
Y u and Yd are 3x3 matrices (i,j are the generation indices) and H is the 
SU(2) doublet Higgs field. After the electroweak symmetry breaking, these 
Yukawa interactions lead to the following quark mass terms 

- L m = ULi{M u )ijU R j + d u (M d ) i:j d Rj + h.c, (2) 
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(M u )ij = (Y v )ijV, 
(M d )ij = (Y d )ijv, 

where v is vacuum expectation value of the neutral component of the 
Higgs field H. 

These mass matrices M u and M d can be diagonalized using two unitary 
matrices U and V 

M u = U u M°Vl (3) 
M d = U d M°Vl (4) 

The diagonal matrix elements (of and M d ) correspond to the exper- 
imentally observed mass eigenvalues. The matrices U and V describe mixing 
between states with definite flavor u, d and with definite mass u', d! for left 
handed and right handed components 

u L = U u u' L , 

u R = V u u' R , (5) 
d L = U d d' L , 
d R = V d d' R . 

The generation mixing in the charged weak current after expressing in 
terms of fields with definite mass is described by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix which consists of two unitary matrices 

Vckm = UtU d . (6) 

In the similar way, assuming that neutrinos are traditional Dirac particles, 
we have mixing for the charged weak current in the lepton sector described 
by Maki-Nakagawa-Sakata (MNS) matrix 

Vmns = U}U V , (7) 
with the relations for neutrino mass M u and charged leptons Mi 

M v = U u M°Vj, (8) 
M [ = UiMfY}. (9) 
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As before diagonal matrix elements of MjP and Mf 1 correspond to the ex- 
perimentally observed mass eigenvalues. 

The present experimental values of the moduli of the quark and leptonic 
mixing matrices are for example given by |2J 



W< 



CKM\ 



I 0.9739 4 0.9751 0.2214 0.227 0.0029 4-0.0045 

0.2214 0.227 0.9730 4 0.9744 0.039 4 0.044 | (10) 
V 0.0048 4 0.014 0.037 4 0.043 0.9990 4 0.9992 



and [H] (similar results given also by [3]) 

/ 0.79 4 0.88 0.47 4 0.61 < 0.20 
|VWs|= 0.19 4 0.52 0.42 4 0.73 0.58 4 0.82 |. (11) 
V 0.20 4 0.53 0.44 4 0.74 0.56 4 0.81 

In some sense matrices M u , M^, M u and Mi are basic quantities because 
they appear in lagrangian after symmetry breaking. If we knew the numerical 
values of their matrix elements we could immediately write the decomposition 
in terms of unitary (ortogonal if there is no CP violation) matrices as in 
eqs.(3), (4), (8) and (9) and then we could calculate Vqkm and Vmns- 

Let us start with some intuitive argument for considering triangular ma- 
trices. One can expect that because the ratios of masses for u-type quarks 
are bigger then for <i-type quarks the mixing for u type quarks is smaller 
then for d type, in other words, that the U u matrix is closer to identity ma- 
trix then Ua- That means in this case that « Vckm- If we also assume 
that Vd ~ I (i.e. Vd is not very different from unity matrix) in fact we look 
in this way for asymmetrical decomposition for Md and we will get in this 
crude approximation from eq.(3) that Md ~ VckmM®. Substituting exper- 
imental values we will get that Md is very close to triangular matrix with 
matrix elements above main diagonal (with very small compared with the 
other matrix elements terms below main diagonal). It seems that in this 
case triangular matrices could be a good starting point to describe quark 
mass matrices. On the other hand we would like to treat quarks and lep- 
tons in the unified way. We assume that triangular matrices describe mass 
matrices not only for quarks but also for neutrinos and charged leptons. In 
this way we have 3 zero elements in the mass matrix, so in some sense we 
follow the papers in which one postulates several zeros in the mass matrix 
jH]. However, we do that in the universal way for quarks and leptons. The 
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advantage is that lower or upper triangular matrices with hierarchical diago- 
nal elements can give very asymmetric decomposition namely that mixing in 
U and V matrices is very different. When we assume matrices in triangular 
form which describe mass matrices of quarks and leptons we have lower and 
upper triangular matrices at our disposal. The second assumption is that 
mass matrices for the members of weak isospin doublets have different tri- 
angular form. That means that when for example d quark mass matrix is 
upper triangular matrix (it has non-zero matrix elements above main diag- 
onal) than u quark mass matrix is lower triangular matrix (non-zero matrix 
elements below main diagonal and zero in the upper part). This assumption 
is connected with the fact that to produce mass of d and u (lower and upper 
weak isodoublet components) with the Higgs mechanism we have coupling 
with Higgs doublet or its hermitian conjugate. The third assumption is con- 
nected with the question how to take into account, treating mass matrices as 
the basic objects, that the right handed components are non active in weak 
interactions. The masses of particles in the weak isomultiplets are very dif- 
ferent. We make an assumption that V matrices responsible for the mixing 
of right handed components in the same weak isodoublet are equal. That of 
course gives a very strong limitation on the possible mass matrices. With 
these three assumptions we are able to calculate mass matrices. Let us start 
with quarks. Instead to solve eq.(3) and (4) for quarks we will numerically 
minimalize the function 



o v {{UlU*)v - Vc X K P Mij ) 2 v (Vuij - V d nf . y (mg - m^f 



r 1 ( m ui m ui ) /in\ 



There are 6 free parameters (3 diagonal and 3 non-diagonal) of upper 
triangular mass matrix and 6 parameters of lower triangular mass matrix 
M u . The values rrf^ = (m d , m s , m b ), ra^ = (m u , m c , m t ) and corresponding 
errors Am^ p , Am„f are taken from the Particle Data Group |2j. When we 
take Vqkm from eq.(10) we would not have in calculations unitarity of Vqkm 
with enough precision and that would result in minimal value of x 2 , so we 
will take elements of Vcku and later V^ns expressed in terms of angles [TU] 
to have unitarity satisfied to high degree and then use calculated in this way 
matrix elements. For the elements of mixing matrices V we take the artificial 
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error 0.01. 

The minimalization of % 2 gives the following results (element of mass 
matrices in MeV) 



M d = 



f 6.16 23.56 14.82 \ 

102.41 177.90 
\ 4246.25 j 

0.9744 0.2250 0.0035 
-0.2249 0.9735 0.0419 
0.0060 -0.0416 0.9991 



(13) 



(14) 



Vl 



0.9999 
0.0132 
5.1 x 10~ 6 




-0.0132 
0.9999 
0.0010 



8.5 x 10~ 6 
-0.0010 

1 



(15) 



(16) 



and 



U u = 






1249.89 
183.07 



1 2.9 x 10~ 5 

-2.9 x 10~ 5 1 
1.4 x 10~ 10 -7.2 x 10~ 6 




M„ 



D 






178000 J 



vl = V] 



(17) 



(18) 



(19) 



(20) 



We see that diagonal matrix elements in mass matrices are very close 
to the "experimental" quark masses for u and d type quarks. Mixing in V d 
is much smaller then in U d and for u quark mixing in U u is smaller then 
in V u . The non-diagonal matrix elements in mass matrices M u and M d are 
comparable what is a little bit surprising. With our assumptions we have got 
the solution for u and d-type quarks in which U u is with high accuracy equal 
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to identity matrix. The Vqkm matrix calculated from eq.(6) is identical to 
Ud obtained in eq.(14). 

We will repeat the calculations for neutrinos and charged leptons. We do 
not know the neutrino masses so it is not possible to solve eqs.(8) and (9). 
To get approximate solutions of eqs.(8) and (9) we will proceed as before and 
minimalize suitably modified y 2 function. In the first two terms in eq.(12) 
quark mixing matrices are substituted by lepton mixing matrices, Vckm by 
Vmns and so on. In the expression for \ 2 given in eq.(12) instead of the third 
term that corresponds to diagonal d quark masses we will take 



Instead of the forth term in eq.(12) for u quark we take the corresponding 
term for charged leptons. The experimental values (Am 2 . 1 ) exp , (Am'^ 1 ) exp 
and corresponding errors are taken from the fit in j^J. Neutrino masses are 
measured in eV while all other masses in MeV. In the case of neutrinos 
equations corresponding to eqs.(13) to (20) are 



(K) 2 -K) 2 -(Am^)^) 2 
(A(Am| 1 ) e ^)) 2 



+ 



((mg) 2 - (mfj 2 - (Am; 
(A(Am| 1 )^)) 2 



2 \exp\2 
31 ) ) 



(21) 




(22) 




(23) 



D 



f 0.00203 





(24) 




/ 0.9901 
0.1403 
V 2.7 x 10~ 5 



-0.1383 
0.9759 
0.169 




(25) 



and 



Mi 



I 0.516 
14.828 106.134 
V -0.096 299.318 1751.51 



(26) 
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u. 



1 

-0.0007 
6.9 x 10~ 6 



MP = 



( 0.511 





0.0007 
0.9999 
-0.0101 



105.658 




17 x 10~ 9 
0.0101 
0.9999 








\ 



1776.99 J 



v} = vl 



(27) 



(28) 



(29) 



We get that (non-diagonal) elements in mixing matrix U u are bigger then 
in V v and in \J\ are smaller then in V\. U\ is not very different from the unit 
matrix. For the Vmns mixing matrix calculated according to eq.(7) we get: 



Vmns 



0.8396 
-0.3800 
0.3881 



0.5432 
0.5875 0.7245 
-0.6000 0.6997 



(30) 



The differences between Vmns and U u matrix elements are smaller then ex- 
perimental errors for matrix elements of Vmns- Non-diagonal matrix ele- 
ments in mass matrix for charged leptons are comparable with those for u 
and d quarks. This is not a case for neutrino mass matrix. Matrix elements 
in V v and in V are not small and comparable with Cabibbo mixing in d- 
type quarks. This result coming from the fit and being a consequence of 
our third assumption is very promising when we think about some ratios 
of charged lepton masses as being comparable in magnitude to the ratios 
of down quark masses. In the neutrino mass matrix non-diagonal matrix 
elements are comparable with the diagonal ones what produces strong mix- 
ing in U v and relatively strong in V v . Diagonal matrix elements in neutrino 
mass matrix eq.(22) are different from the diagonal masses in the matrix 
eq.(24) because of strong mixing. According to eq.(24) we get: 



mi„ = 0.002 eV, 

m 2v = 0.009 eV, (31) 
m 3u = 0.048 eV. 

These values of neutrino masses come out from the fit. We get some 
hierarchy in neutrino mases but not comparable with what we have for quarks 
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and charged leptons. Specially interesting is the first mass eigenvalue because 
it sets a scale. Its mass is about 1/4 of the mass of the second neutrino. The 
values of m,2u and m^ v are close to what one would expect from mass squared 
differences with small first neutrino mass. The question is what is the error 
for m\ v in the fit. When we fix m\ v and make fit for other quantities we get 
minimal x 2 smaller then one for 0.002 < m,\ v < 0.005. For much smaller 
values of m\ v we also get Xmin < 1 so we do not have strong limitations 
for m\ v . It seems that except of neutrino when mixing in mass matrix is 
very strong and the scale of masses is completely different the non-diagonal 
matrix elements in mass matrices of u and d type quarks and charged leptons 
are comparable (not very different in magnitude) in spite of the differences 
in ratios of masses of quarks and charged leptons. Could it be a hint in 
understanding the structure of masses of quarks and leptons? 

Till now all the mixing matrices were real (orthogonal). There is no 
problem with repeating the calculation taking into account CP violation. In 
this case we simply have additional imaginary parts in Vckm matrix and 
unitary instead orthogonal mixing matrices. When we repeat the calculation 
as before for u and d type quarks solving eqs.(3) and (4) by minimizing x 2 
function corresponding to that given in eq.(12) and taking into account that 
we have complex quantities we get: 




6.16 23.55 - O.Oli 
102.41 





(32) 



0.9744 0.2250 
-0.2249 - O.OOOli 0.9735 - 3 x 10-H 
0.0077 - 0.0029? -0.0412 - 0.0007i 



0.0018 -0.0030i \ 
0.0419 

0.9991 j 



U d = 



(33) 



V 




D 



( 



6 
105 
4250 



(34) 




/ 0.9999 -0.0132 

0.0132 0.9999 
V (0.3 + OAi) x 10~ 5 0.0011 + 1.9 x 10 5 i 



(1.1 + 0.4i) x 10~ 5 
0.0011 + 1.7 x 10- 5 i 
1 
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and 




( 2.75 

16.49 - 0.0002^ 1249.89 j (36) 

V 0.47 + 0.77^ 181.41 + 2.96z 



(4.1 - 6.7i) x 10~ n 
(7.8-O.lz) x 10~ 6 
1 

(37) 



/ 


1 


2.9 x 10 


-5 




-2.9 x lO" 5 


1 




V (1-6 


-0.6i) x 10~ 10 


(-7.2-0.li) 


x 10~ 6 






/ 2.75 









1250 









v 


178000 






V* = vl 

u a 





(38) 



(39) 

Eqs.(32) to (39) should be compared with those given in eqs.( 13-20). 
We see that there is no big change except of additional imaginary parts. 
There is no problem with repeating this type of calculation for neutrinos and 
charged leptons taking into account CP violation but we do not have enough 
information about Vmns matrix at the moment. 

We have shown that it is possible to determine mass matrices for quarks 
and leptons under the assumption that they are triangular and that the 
matrices for upper and lower weak isospin components have zeros below or 
above main diagonal. The mixing matrices corresponding to diagonalization 
of weakly non-interacting right handed components are assumed in weak 
isodoublet to be equal. As an input the values of matrix elements for CKM 
and MNS matrices with known errors and masses of quarks and charged 
leptons and in the case of neutrinos squared mass differences also with corre- 
sponding errors are used. The masses of neutrinos are determined in the fit 
and the possible range of smallest neutrino mass is given. The U matrices for 
up type quarks and charged leptons are not very different from unit matrix. 
The V matrix for down type quarks is close to unit matrix and for neutrino 
matrix elements of mixing matrix are not small and comparable with Cabibbo 
mixing. In this solution U matrices for d type quarks and neutrinos are very 
close to measured Vckm and Umns matrices. It was actually our aim to find 
out whether it is possible to get this type of solution. For up, down and 
charged leptons the non-diagonal matrix elements in the mass matrices are 
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of the same order what is not expected taking into account masses of these 
particles (up and down type quarks and charged leptons). Neutrino masses 
have different scale (the hierarchy is rather weak) and some non-diagonal 
matrix elements are comparable with the diagonal. In the case of up and 
down type quarks we also calculated using our method mass matrices and 
mixing matrices taking into account CP violating phase. 
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